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USE  OF  TIME-DEPENDENT  CHEMICAL  SENSOR  SIGNALS 
FOR  SELECTIVE  IDENTIFICATION 


1.  INTRODUCTION 

The  magnitude  of  a  sensor  response  at  steady  state  is 
commonly  used  to  determine  the  species  concentration;  whereas, 
the  information  in  the  transient  response  is  generally  not 
processed.  However,  the  transient  signal  can  be  made  to  contain 
information  that  can  identify  and  quantify  the  chemical  species 
present.  In  this  report,  a  purposely  caused,  periodic, 
nonsteady-state  sensor  signal  is  generated,  and  this  modulated 
signal  is  controlled  experimentally  to  provide  information  about 
specific  chemical  species  present  and  their  concent  rat  ions . i  i 

Modulated  systems  have  been  employed  in  multiplex  gas 
chromatography  in  which  multiple  samples  are  injected,  the 
column  temperature  is  modulated,  and  the  cross-correlation  of 
the  resulting  chromatogram  is  computed. <-6  a  variety  of 
modulated  detectors  have  been  reported  recently,  including  a 
mass-flow  modulated,  flame  ionization  detector, a  temperature- 
modulated,  thermal  conductivity  detector;^  a  temperature- 
modulated,  array  of  MOS  gas  detectors;^  and  a  solute-modulated, 
electron  capture  detector. lo  jn  all  of  these  detectors,  the 
information  content  of  the  complex  signal  is  higher,  and  the 
speed  of  data  acquisition  is  greater  than  in  corresponding 
signals  obtained  using  conventional  systems. Synchronous 
detection  of  modulated  signals  also  can  improve  the  signal-to- 
noise  (S/N)  ratio  because  low  frequency  noise  and  drift  are  not 
modulated  and,  therefore,  can  be  eliminated  from  the  amplified 
signal . 


A  theory  has  been  developed  for  an  existing  sensor 
system, 1"3  and  measurements  were  made  to  aid  in  model  develop¬ 
ment.  The  system  consists  of  a  reaction  chamber  and  an 
electrochemical  ceil  (Figure  1).  Tne  gas  to  be  analyzed  i:^ 
exposed  to  the  electrically  heated  catalytic  filament  in  the 
reaction  chamber  that  precedes  the  detectoi  .  The  rate  at  whic'n 
the  chemical  reaction  proceeds  depends  on  the  species  concentra¬ 
tion,  catalyst  activity,  and  chamber  residence  time.  After 
leaving  the  reaction  chamber,  the  products  of  reaction  are 
detected  by  an  electrochemical  sensor  tha 
proportional  to  the  product  concentration.  The  product 
concentration  can  be  controlled  most  easily,  in  this  c'ase,  by 
changing  the  temperature  of  catalyst  or  the  flow  tale  (i.e., 
residence  time  of  the  analyte  in  the  reactor). 


7 


Potwnioalal 


Ga* 

Air 


Figure  1.  Detector  System. 


The  transient  response  of  the  sensor  is  expressed  in 
terms  of  the  volume  (Vq)  of  the  reaction  chamber,  the  reaction 
rate  constant  (k),  the  concentration  [C(t)]  of  the  molecules 
detected  by  the  sensor,  the  time  constant  of  the  electrochemical 
sensor  (is) f  and  the  absolute  sensor  sensitivity  (m)  in  the 
microamperes  per  parts  per  million.  Measurements  were  made 
under  two  different  conditions  using  low  concentrations  of 
cyclohexane  and  benzene  in  air.  One  condition  is  with  a  fixed 
filament  temperature,  and  the  other  condition  is  with  a  periodic 
modulation  of  the  filament  temperature,  producing  a  constant 
concentration  and  a  periodic  modulation  in  the  concentration  of 
the  reactant  species,  respectively.  By  making  transient 
measurements  at  different  flow  rates,  the  parameters  character¬ 
izing  the  system  response  were  determined,  and  models  for  the 
system  operation  were  developed. 


2.  THEORETICAL  ANALYSIS 

2 . 1  Model  of  Electrochemical  Sensor. 

An  equivalent  circuit  model  proposed  for  the  electro¬ 
chemical  sensor  is  shown  in  Figure  2.  This  model  assumes  that 
the  electrochemical  cell  is  a  linear  amperometric  sensor  with  a 
steady-state  output  current  (is(t)]  proportional  to  the 
concentration  [C(t)]: 


i,,(0  =  "» 
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and  tne  current  [i(t)l  measured  by  the  potentiostat  detecLoi 
circuit  is  proportional  to  the  voltage  [V(t)]: 


]  ^  V{t)xG 


(  ib) 


where  G  is  constant.  Implicit  in  equation  lb  is  the  assuirpiion 
that  the  response  time  of  the  potentiostat  is  much  shorter  than 
the  response  time  of  the  electrochemical  cell.  The  genetal 
solution  for  i(t)  is 


i{t)  = 


r, 


(  2  ) 


where 

C(t)  =  0  for  t  <  0 

m  =  m'RG  =  intrinsic  sensitivity  of  cell 
Is  =  RC  =  sensor  time  constant 


If  the  number  of  moles  of  gas  in  Vq  increases  from  0 
t  =  0,  the  resulting  current  response  is 

i(t)  ■-  m(n,/ 


to  nc  at 


(  3a) 


WE 


Figure  2.  Equivalent  Circuit  for  the  Elect  r  ochem  i  ca  1  t’c  :  i  . 

CE  is  the  counter  electrode;  WE  is  the  woikii.g 
electrode;  is(t)  is  the  Faradaic  current;  K  and  C 
are  the  internal  resistance  and  capacitance  of  the 
elect rochemical  cell. 


One  interesting  consequence  of  this  circuit  modc'l  is 
that  the  potentiostat  output  i(t)  is  directly  pr '  ,p. ' t  t  i  ona t> 
the  instantaneous  concentration  C(t)  whenever  c.  ;(t  )  ct  h  . 
This  condition  occurs  at  steady  state  as  well  ac  at  t  r.L  ^arpat 
peaks  and  valleys  that  occur  during  periodic  no  du  1  at  i  a  n.  d  tee 
analyte  concentration.  To  prove  that  the  output  is  diiHctly 
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proportional  to  the  concentration,  we  note  that  the  current 
is{t)  in  the  circuit  shown  in  Figure  1  is  given  by 


«,,(n  =  — H — X - 

^  ^  R  C  dt 


(3b) 


When  di(t)/dt  =  0,  it  follows  from  equation  lb  that  dV ( t ) / 
dt  -  0.  Therefore,  no  current  flows  through  the  capacitor,  and 
from  equation  3b,  it  follows  that 


*40  = 


R 


(3c) 


Using  this  result  and  equation  la,  we  can  find  that  the  sensor 
output  will  be  pr  .portional  to  the  concentration  whenever 
di ( t ) /dt  =  0 : 


i[t)  =  mC{t)  when  di[t)l dt—Q 


(3d) 


This  condition  occurs  at  steady  state  or  at  peaks  and  valleys  in 
the  modulated  sensor  signal. 

Figure  3  shows  the  effect  of  ts  on  i(t)  when  the  input 
concentration  is  varied  sinusoidally  at  a  frequency  u .  The  peak 
at  which  output  current  shifts  from  the  peak  of  the  input  by  an 
amount  equal  to  tan-i(cots),  and  the  output  current  is  reduced  by 
the  factor  [1  +  ( cots ) 2  ] -i/2 .  Curve  (a)  in  Figure  3  tracks  the 
input  concentration  almost  exactly.  This  model  predicts  that 
the  sensor  response  time  must  be  significantly  smaller  than  the 
period  of  modulation  to  insure  an  adequate  signal.  Faster 
sensors  will  allow  higher  modulation  frequencies  for  given 
signal  amplitudes.  The  electrochemical  sensor  may  not  be  an 
ideal  detector  in  their  kind  of  sensor  system. 

2.2  Model  of  Reaction  Chamber  Dynamics. 

Input  gas,  with  component  species  A  and  B,  flows  into 
the  reaction  chamber  of  Vq  at  a  rate  of  u  (liters  per  second). 

A  reaction  takes  place  in  which  a  reactant  C  is  formed 


k 

A  +  B  -<•  C 

where  k  is  the  reaction  rate  constant  (liters  per  mole-second). 
The  gas  in  the  chamber,  which  is  assumed  to  be  homogeneous  and 
isothermal,  exits  into  an  electrochemical  sensor  where  it  is 
detected.  The  chamber  contains  a  fixed  amount  of  catalyst .  it 
is  assumed  that  the  gas  has  returned  to  the  ambient  temperature 
when  it  reaches  the  sensor  and  that  the  flow  through  the  sensor 


Timt  Cytt, 

Curvt  P«rio4  PtokArtipl. 

a  0.01  .t»« 

b  0.1  JIT 


Figure  3.  Elect rochemical  Cell  Output  for  a  Periodic  Variation 
in  Concentration  with  Frequency  to  =  2ii/T  for 
Different  Ratios  of  the  Cell  Time  Constant  to  the 
Period  is/T:  (a)  0.01,  (b)  0.1,  (c)  0.3,  and 
(d)  1.0. 


The  reaction  rate  r(t)  (moles  per  liter-second)  is 
directly  proportional  to  the  concentrations  of  the  reactants: 


r  = 


(4) 


In  the  typical  trace  pollutant  detection  problem,  the  concentra¬ 
tion  of  the  species  to  be  detected  (species  A)  is  several 
hundred  parts  per  million  whereas  the  oxygen  concentration  is 
200,000  ppm.  To  solve  the  system  of  coupled  first-order, 
nonlinear  differential  equations  describing  the  process  in  the 
reaction  chamber,  we  assume  the  following: 

•  The  change  in  the  concentration  of  oxygen  due  to  the 
reaction  is  negligible. 

•  The  nitrogen  is  an  inert  diluent  gas  and  does  not 
affect  the  form  of  the  equations. 

•  The  flow  of  air  into  Vq  equals  the  flow  out  (total 
flow  in  =  total  flow  out). 

The  expressions  obtained  for  iic  ( t )  ,  the  number  of  moles 
of  species  C  in  Vq ,  are  parameterized  in  terms  of  the  residence 
time  Xres  and  the  characteristic  chemical  reaction  time  ta  where 


1 1 


u 


(  5a  ) 


and  u  =  flow  rate  of  sample  gas  into  chamber  (liter  per  second) 
and 


Ta  =  Vol{kn  b) 


The  reaction  rate  (r)  can  be  written  in  terms  of  Xa-' 


r  =  («a/^o)(lAa) 


The  reaction  rate  constant  (k)  is  assumed  to  depend  on  tempera¬ 
ture  only: 


k  =  A  exp{-AE/R  ^T) 


where 


A  =  pre-exponential  factor  for  the  particular  oxidation 
AE  =  activation  energy  for  conversion  of  A  to  C 
Rg  =  gas  constant  per  mole 

In  the  Arrhenius  equation,  the  pre-exponential  typically  has  a 
square  root  dependence  on  temperaturei2  (neglected  in  our 
analysis ) . 

In  the  following  paragraphs,  solutions  for  nc(t),  the 
modulated  concentration,  are  determined  for  two  cases. 

2,2.1  Case  1  (Constant  Temperature), 

When  the  filament  in  the  reaction  chamber  is  held  at  a 
constant  temperature,  for  a  step  function  change  in  the  number 
of  moles  of  species  A  into  Vg  at  time  t  =  0 ,  we  find: 


n,{t)  =  ) 


where 


(t-*)  ^  =  {rrf>]  ‘  + 


The  equation  for  (x*)-i  can  be  interpreted  as  representing  the 
sum  of  the  number  of  moles  per  second  due  to  the  flow  plus  the 


number  of  moles  per  second  due  to  the  reaction.  The  steady- 
state  solution  (as  t  -►infinity)  is: 


njoo) 


(oo)  =  k 


(10a) 

(10b) 


2.2.2  Case  2  (Time-Dependent  Filament  Temperature). 

Variation  of  the  temperature  with  time  produces 
exponential  variations  in  the  reaction  rate.  The  following 
general  integral  expression  for  n(t)  can  be  written  for  this 
case : 


n(0 


f r(t)dt 


(11) 


2 . 3  Model  for  the  Detection  System. 

The  model  for  the  detection  system  is  developed  througli 
a  joint  consideration  of  the  sensor  model  response  and  the  time- 
dependent  production  of  analyte  C  in  the  reaction  chamber.  The 
expression  for  the  output  current  from  the  electrochemical 
sensor  is  obtained  by  substituting  the  appropriate  equation  for 
the  concentration  of  the  oxidation  product  C  into  equation  1  for 
the  sensor  response.  For  the  constant  temperature  case,  with  a 
step  function  change  in  the  number  of  moles  of  A  in  the  reaction 
chamber  at  t  =  0,  we  obtain: 


■(')  = 

*  I'd  *  o 


(12) 


The  first  term  in  parentheses  on  the  right-hand  side  of  the 
equation  represents  the  fraction  of  A  converted  to  C  [the 
conversion  ratio  (R)].  The  second  term  represents  the  initial 
concentration  of  A  in  the  reaction  chamber,  and  the  last  term 
yields  the  time  dependence.  At  steady  state,  this  last  term 
goes  to  unity  as  t  -►  The  conversion  ratio  (R)  can  be  written 

using  equation  5: 


/?  -  (1  f  u/(A:n  i,))  ‘ 


(13) 
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Note  that  R  depends  on  the  flow  rate  (u)  and  the  reaction  rate 
(k)  but  is  independent  of  cell  volume  (Vq).  Therefore,  steady- 
state  measurements  can  give  information  about  m  and  k  but  not 
about  ts  or  Vq.  The  constant  m  is  characteristic  of  the  sensor 
sensitivity  for  a  given  oxidation  product,  and  k  is  characteris¬ 
tic  of  the  analyte  for  a  given  catalyst  and  reaction  chamber 
condition. 

To  obtain  values  for  these  parameters,  we  fit  solutions 
for  i(t)  to  measurements  taken  at  various  flow  rates.  At  very 
low  flow  rates,  R  goes  to  unity,  and  all  of  the  reactant  A  is 
converted  to  C.  Conversely,  at  high  flow  rates,  the  conversion 
is  low,  and  R  is  given  by: 


R  =  (fcn(,)(l/u),  for  u  large 


(14) 


Therefore,  plots  of  the  steady-state  current  i(‘»)  can  be  used  to 
identify  the  region  of  operation  and  the  behavior  of  the  system 
as  shown  in  Figure  4.  The  constants  m  and  k  can  be  obtained 
from  the  steady-state  data  by  plotting  the  inverse  of  the 
steady-state  signal  per  unit  concentration  (  na/Vo)  [  1/i  («>)  ] 
versus  u  since 


(»^/j/'^<.)(l/»(oo))  =  h)  (15) 


Data  acquired  at  a  low  flow  rate  can  be  used  to  determine  m.  To 
determine  the  activation  energy  (AE),  measurements  must  be  made 
at  two  or  more  different  temperatures.  The  remaining  parameters 
Vo  and  Is  must  be  determined  by  fitting  transient  measurements 
to  the  theoretical  expressions  for  i(t).  The  mathematical 
expressions  are  simplified  if  the  measurements  are  made  in 
either  a  fast  or  slow  flow  regime. 


Figure  4.  Steady-State  Output  Current  Versus  Inverse  of  the 
Flow  Rate  u. 
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2.4 


Use  of  Periodic  Filament  Temperature  Variations  in  the 
Detection  System. 

From  the  above  consideration,  we  are  now  prepared  to 
investigate  how  a  periodic  variation  (e.g.,  sinusoidal)  in  the 
filament  temperature  can  yield  information  that  can  assist  in 
differentiating  between  different  species  that  produce 
detectable  reaction  products  (C)  from  reactions  with  differnnc 
rates  (e.g.,  different  rate  constants  and  activation  energies). 
For  simplicity  in  this  analysis,  we  assume  that  the  variations 
in  temperature  T  about  Tq  that  are  given,  for  example,  by 

T  =  AT  sin(yjt+4>)  T„  (lb) 


are  sufficiently  slow  so  that  the  solution  for  the  steady-state 
current  (equations  15  and  13)  is 


j’(cx))  =  mC^R 


(17) 


where  Ca  =  na/Vo  is  valid  for  all  times  during  the  temperature 
cycle.  This  condition  will  be  met  if  the  period  of  the 
modulation  is  much  greater  than  i* . 

We  can  write  a  two-term  Taylor  series  expansion  for  the 
steady-state  current  as  a  function  of  the  filament  temperatiue 
(T) 


.(oo,r)=,(oo)|r. +  ^!^lr„Ar 


(18) 


The  first  term  on  the  right  represents  the  dc  current  (idc)- 
The  second  term  represents  the  ac  current  magnitude  (iac)  in  the 
small  signal  approximation  in  which  all  higher  order  terms  in  At 
have  been  omitted  in  equation  18.  This  approximation  will  be 
valid  provided 


1 


di (oo) 


AT  «  I 


hie  »(oo)  dT 

Differentiating  equation  17,  we  obtain 


{  19  ) 


AT 


di{oo) 

dT 


T„ 


mCnR  (1  R) 


AE  AT 
R,;  T  T 


(  20  ) 
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Because  0  <  R  <  1  for  any  finite  nonzero  flow,  it  follows  that 
the  derivative  in  equation  20  never  vanishes;  thus,  there  is 
always  an  ac  response.  Because  i(®)  =  idc/  from  equation  18,  it 
follows  that 


7^  =  (i-fi) 

*(ic 


AE  AT 
RgT  T 


(21) 


This  ratio  is  independent  of  concentration  and,  therefore, 
carries  the  qualitative  information  required  to  identify  the 
reacting  species  A.  The  ratio  depends  only  on  the  flow  rate. 
At,  and  AE.  Therefore,  it  should  be  possible  to  distinguish 
between  reactions  and  components  with  different  activation 
energies . 


The  magnitude  of  the  dc  signal  varies  directly  with  R 
(equation  17);  whereas,  the  ratio  of  the  ac  to  dc  signal  goes  as 
l-R.  Therefore,  at  low  conversions,  the  dc  signal  is  reduced, 
but  the  ratio  of  the  ac  to  dc  signal  is  increased.  For  a 
typical  case  (AE  =  30  kcal/mole,  T  =  900  K,  T  =  10  “C),  we  find 
that  iac/idc  varies  from  0.18  for  10%  conversion  to  0.10  for  50% 
conversion . 


The  operating  point  [Pq  (u,T)]  at  which  the  ac  signal  is 
a  maximum  can  be  determined  by  computing  a  second  derivative  and 
setting  it  equal  to  zero 


dT 


(i^i(oo) 

dT'^ 


P  =0 

*  o 


(22) 


This  yields  a  condition  on  the  dc  signal  for  which  the  ac  signal 
is  a  maximum: 


«{oo)  |p„ 


mC„  ^  2R(]  T 

^  AE 


(23) 


The  peak  in  the  ac  curve  corresponds  to  a  dc  conversion  ratio 
slightly  less  than  50%: 


R 


P,, 


1  RgT 

2  AE 


(24) 


For  the  previous  numerical  example  R 1 Po  =  0.45.  The  physical 
significance  of  equation  24  becomes  apparent  by  noting  from 
equations  5a  and  5b  that 


Ta 


T 


res 


U 

knt, 


(25) 
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Therefore,  at  Pq/  it  follows  that 


r„  1/2+RgTI/^E 

-±-  ^  -i - -  >1  (26^ 

r,,,  II2~RgTI/^E  - 

For  most  practical  values  of  T  and  AE,  it  follows  that  la  “  ires 
at  Pq.  Therefore,  the  ac  signal  is  at  its  peak  when  the 
characteristic  reaction  time  approximately  equals  the  residence 
time . 

2 . 5  Analysis  of  Gas  Mixtures. 

If  there  are  independent  noninterferr ing  species  present 
and  the  system  response  to  each  is  linear,  gas  mixtures  can  be 
analyzed  by  the  described  detector  system.  It  is  more  likely 
that  such  approximations  will  be  valid  if  low  concentrations  and 
conversion  ratios  are  maintained. 


2.5.1  Constant  Filament  Temperature. 

A  sensitivity  Sij  (ampere  per  parts  per  million)  for  the 
species  i  at  temperature  Tj  and  flow  rate  uj  can  be  defined  as: 


c  _  I 

I  T,,u, 


(27) 


where  Ci  =  ni/Vo  is  the  concentration  of  species  i.  General¬ 
izing  our  previous  result  (equation  15)  yields: 


(  28) 


If  the  terra  Sij  is  determined  for  each  species,  under  a  set  of 
flow  rate  and  temperature  conditions,  the  total  steady-state 
response  at  condition  j  is 


Ij[oo)  —  Yi  ^t^i]  i  21i ) 


To  solve  for  the  concentration  of  each  species  {C^,  i  -  l,...n), 
calibration  measurements  must  be  made  to  determine  Sj,  ar  n 
different  sets  of  conditions  j  =  l,...n.  The  sensitiviry  Sj,  is 
a  square  matrix  that  can  be  inverted  to  find  the  Ci-  To  invert 
S,  the  coefficients  for  Sij  must  be  different  for  different 
species  (i  =  l,...,n)  at  the  same  conditions  (fixed  j)  and  be 
different  for  the  same  species  (fixed  i)  at  different  condi¬ 
tions.  Also,  elements  of  Sij  must  be  large  enough  to  permit 
accurate  computations. 
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2.5.2  Periodic  Variation  of  Filament  Temperature. 

Generalizing  equation  20  to  several  independent  species 
indicated  by  the  subscript  “i,"  and  using  equation  13  for  R,  we 
obtain  an  expression  for  the  total  ac  response  at  conditions 
{ u j  /  Tj )  : 


»«C,M  =  Arsc.s.y 

i 


(30) 


where 


AE, 

ulk^Tib 

Ra 

(u/kini,  +  1)2 

(31) 


The  factor  sij  is  dependent  on  the  temperature,  u,  and  AE^.  To 
ensure  that  we  can  solve  the  system  of  equations,  Sij  must  obey 
the  same  restrictions  as  Sij. 


3.  MODEL  VALIDATION 

Figure  5  shows  the  response  of  the  electrochemical 
sensor  to  20  ppm  CO/air  gas  mixture  (CO  sensor,  ND,  model  110, 
Pittsburg,  PA).  The  sensor  gave  a  response  time  constant  of 
5.71  s,  a  sensitivity  of  3.75  A/ppm,  and  operates  consistently 
with  the  model  proposed  (Figure  1  and  equation  3a).  Data  for  CO 
taken  at  different  flow  rates  agreed  with  equation  15.i'3  In 
Figure  6,  the  response  of  the  electrochemical  sensor  is  shown  at 
different  frequencies  of  the  filament  temperature  cycles  of  the 
reactor  filament  when  a  benzene  air  mixture  is  introduced.  The 
amplitude  of  the  temperature  cycle  is  from  650  to  870  °C.  The 
oscillation  is  damped  due  to  the  slow  sensor  response  (Figure  3) 
or  to  the  slow  filament  response  at  the  higher  frequencies.  At 
low  frequencies,  the  steady-state  response  i('^)  is  nearly 
achieved  at  each  temperature  extreme. 


Figure  5.  Sensor  Response  to  CO.  Solid  line  is  best 
exponential  fit  of  the  form  of  equation  3a. 
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An  experiment  similar  to  that  shown  in  Figure  6  was 
performed  for  HCN  in  which  the  amplitudes  at  the  top  and  bottom 
of  the  peaks  (called  "A"  and  "B,"  respectively)  were  studied.  A 
consequence  of  the  sensor  model  is  that  the  signal  should  be 
proportional  to  the  concentration  whenever  di/dt  =  0  (Section 
2.1).  This  condition  is  met  at  the  peaks  and  valleys.  The 
model  is  validated  by  the  observations  (Figure  7),  and  linear 
response  is  obtained  for  this  sensor  system  even  though  it  is 
not  at  steady  state. 


FReOUeNCY(Hi) 

0.3  0.4  0  3  0  2  0.1  0  03  0  025 


Figure  6.  Electrochemical  Cell  Response  to  511  ppm  Benzene/Air 
at  Different  Filament  Frequencies  with  Square  Wave 
Power  Input,  Temperature  Variation  Between  650  and 
870  '’C,  u  =  120  cm3/min,  Vq  =  1  cm3. 


Figure  7.  Linearity  of  Transient  Part  of  Sensor  Response  Versus 
Concentration  of  HCN.  "A"  and  "B"  represent  peaks 
and  valleys,  respectively. 
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The  steady-state  detector  response  versus  temperature  is 
shown  in  Figure  8  for  512  ppm  of  benzene.  The  curve  represents 
the  best  fit  of  equation  15  to  the  first  5  data  points.  The 
drop  in  response  at  high  temperatures  is  not  predicted  by  the 
theory,  but  the  theory  is  valid  at  the  lower  temperatures.  A 
reductio  .  in  oxidation  rate  at  increasing  temperature  suggests  a 
more  complex  reaction  mechanism  than  can  be  explained  by 
Langmuir-Hinshelwood  kinetics. 13  The  theory  could  be  extended 
to  include  these  mechanisms,  but  it  is  easier  to  restrict  our 
discussion  to  temperatures  below  700  “C  where  our  model  is 
operative. 


Figure  8. 


Current  Versus  Temperature  Data  for  511  ppm  Benzene, 
The  solid  line  is  the  best  fit  to  equation  15. 


Assuming  the  data  for  benzene  is  representative  of  other 
species,  equation  15  can  be  used  to  simulate  the  steady-state 
sensor  response  to  changes  in  the  reaction  rate  due  to  changes 
in  AE  (i.e.,  species)  and  the  filament  temperature  (i.e.,  rate 
constant).  In  Figure  9,  the  calculated  steady-state  current  is 
plotted  as  a  function  of  temperature  for  three  different 
chemical  species  that  have  three  different  AE.  A  response  of 
50/ijA  represents  100%  conversion  of  A  to  C  (m  =  0.5  A/ppm).  The 
position  of  the  i  versus  T  curves  is  very  dependent  on  Ae.  This 
dependency  suggests  the  possibility  of  discrimination  between 
species  with  different  reaction  rates. 
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u  -  1  cm3/min  and  A  =  1.5  x  106  cm3/niole-s 


Figure  9.  Calculated  Current  Versus  Temperature  Curves  for 
Three  Different  Activation  Energies  (AE),  u  = 

1  cm3/min,  A  =  1.5  x  106  cm3/mole/s. 


Figure  10  illustrates  a  simulation  of  the  sensor 
response  to  a  mixture  of  species  with  the  AE  (AE  =  25,  35,  and 
45  kcal/mole)  used  in  Figure  9.  Because  it  is  assumed  that  the 
chemicals  give  independent  responses,  Figure  10  is  simply  the 
sum  of  the  individual  curves  from  Figure  9.  When  the 
derivative,  with  respect  to  temperature,  is  taken  of  the  mixture 
curve,  a  "spectrum"  is  obtained  with  peaks  corresponding  to 
different  activation  energies  (Figure  11).  The  peak  heights 
correspond  to  the  magnitude  of  the  ac  sensor  response,  and  the 
location  of  the  peak  identifies  Pq  and  To  at  whicn  the  sensor 
exhibits  the  greater  temperature  dependence  for  the  species  with 
that  particular  AEq  (equations  20  and  22).  This  kind  of 
"spectrum"  could  be  generated  experimentally  by  modulating  the 
temperature  in  a  sine  wave  while  also  ramping  the  temperatur^r 
from  low  to  high  and  using  phase  sensitive  detection  to  monitor 
only  the  modulated  sensor  response. 

Modulation  of  the  filament  temperature  about  !'□  should 
not  produce  a  modulated  sensor  signal  for  spvecies  with  AF 
significantly  different  from  AEq.  This  is  illustrated  in  Figure' 
12  where  the  temperature  is  modulated  in  a  sine  wave  about  Tq 
for  a  species  with  a  AEq  of  33.7  kcal/mole  (and  A  -  1.5  x  106 
cm3/mole/s).  The  calculated  sensor  response  shows  very  little 
ac  component  for  compounds  with  AE  of  25  and  40  kcal.  ino'e.  The 
absolute  current  scale  depends  cn  the  sensitivity  of  the  sensir. 
The  conversion  rate  of  AE  =  25  kcal/mole  is  0.99;  and  at  AF  - 
33.7  kcal/mole,  the  ac  signal  is  a  maximum. 


21 


60 


T=700+SIN<PI»l)t*10  CC) 


&E*23 


0  2  5  5  7.5  10 

I  CSEC) 


Figure  12.  Calculated  Response  of  Sensor  for  Activation 

Energies  from  25  to  40  kcal/mole.  The  filament 
temperature  varies  sinusoidally  from  690  to 
710  °C  with  a  period  of  2  s. 


The  simulations  have  assumed  that  the  pre-exponential 
term  is  the  same  for  each  species,  and  the  differences  in  rate 
constants  are  due  solely  to  differences  in  activation  energies. 
This  is  obviously  not  always  a  valid  assumption.  One  can  also 
expect  the  pre-exponential  term  A  to  vary  from  species  to 
species,  which  will  also  cause  the  current  versus  temperature 
curves  to  shift.  However,  the  greatest  discrimination  with 
temperature  variation  will  result  from  differences  in  AE  becdust' 
AE  appears  in  the  exponent  of  the  rate  constant.  The  ac  simula¬ 
tions  assume  that  temperature  variations  are  sufficiently  slow 
that  equation  15  remains  valid  (i.e.,  steady-state  currents  art- 
observed).  This  assumption  was  made  for  mathematical  simplicity 
and  could  be  el  i..unated  by  use  of  equation  11. 


4.  CONCLUSIONS 

A  model  has  been  constructed  to  analyze  the  output 
current  from  a  detection  system  that  includes  a  modulat  u  t 
concentration  and  a  detector.  Specific  responses  hu-.u-  br-  i 
discussed  for  a  system  consisting  of  a  reaction  chambe;  wj i n  a 
heated  catalytic  filament  and  an  electrochemical  sensor.  llie 
parameters  required  to  apply  the  theory  can  be  obtained  t'y 
system  operation  at  different  flow  rates  and  t empe r a t u  lou- . 
Theory  predicts  that  periodic  modulated  reaction  ta.es  can  ue 
use  in  identifying  and  quantifying  gases  alone  and  in  mixtiiifu. 
Initial  data  for  our  detector  system  with  HCN  and  iien.'ene  diul 
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computer  simulations  confirm  certain  aspects  of  the  model  and 
outline  its  limitations.  In  particular,  it  appears  possible  to 
determine  concentrations  using  appropriate  measurements  taken 
before  steady  state  is  reached.  Therefore,  by  using  the 
approach  proposed,  the  response  time  for  concentration  measure¬ 
ments  from  a  given  sensor  can  be  reduced.  A  "smart"  sensor  that 
cycles  rapidly  at  high  concentrations  can  immediately  detect 
dangerous  situations;  whereas,  one  that  cycles  slowly  at  low 
concentrations  achieves  maximum  sensitivity.  The  measurement 
accuracy  may  be  reduced  using  ac  methods  if  the  dc  signal/noise 
ratio  is  excellent.  On  the  other  hand,  if  significant  low 
frequency  noise  or  drift  is  present,  the  use  of  synchronous  ac 
detection  may  be  more  accurate. 

More  extensive  experimentation  will  be  required  to 
validate  and  refine  the  model  so  that  practical  applications  can 
be  realized  for  various  concentration-modulated  systems.  Some 
potential  advantages  of  using  nonsteady-state  responses  from  a 
detector  for  chemical  analysis  include  faster  response,  greater 
dynamic  range,  and  the  ability  to  carry  quantitative  and 
qualitative  information  about  several  species  on  a  wave  form 
type  output.  The  filament/catalyst  and  electrochemical  sensor 
system  may  be  too  slow  to  explore  the  full  potential  of  the 
approach  outlined  herein. 

The  computer  simulations  and  theoretical  analyses  appear 
to  indicate  that  the  appropriate  choice  of  filament  temperature, 
amplitude  of  temperature  modulation,  and  flow  rate  will  result 
in  the  "filtering  out"  of  responses  from  all  species  with 
significantly  different  rate  constants  than  that  of  a  selected 
compound  to  be  detected  by  the  sensor  system.  Experiments  are 
now  required  to  determine  the  actual  variability  of  k  and  AE 
from  species  to  species.  A  more  general  sensor  response  model 
should  be  incorporated  into  the  theory,  and  faster 
modulator/detector  pairs  should  be  studied  to  validate  the  above 
arguments . 
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